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Abstract 

The principal pivot transform (PPT) of a matrix A partitioned relative to an invertible 
leading principal submatrix is a matrix B such that 

A ( Xl ) = ( Vl ) if and only if 



X2 J \ 2/2 / \ X 2 J \V2 

where all vectors are partitioned conformally to A. The purpose of this paper is to survey 
the properties and manifestations of PPTs relative to arbitrary principal submatrices, make 
some new observations, present and possibly motivate further applications of PPTs in matrix 
theory. We pay special attention to PPTs of matrices whose principal minors are positive. 
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1 Introduction 



Suppose that A € M n ((D) (the n-by-n complex matrices) is partitioned in blocks as 

W* 1 A ") (1.D 



,^21 ^22. 

and further suppose that j4h is an invertible submatrix. Consider the matrix 

B=( {All) ~ l -{An)- l A 12 \ 

\A 2l {A ll )- 1 A 22 -A 21 (A 11 )- 1 A 12 J- { - } 

The matrices A and B are related as follows: If x = (xJ,x 2 ) T and y = (yi,y 2 ) T in (D n are 
partitioned conformally to A, then (see Theorem |3.1|) 

A^W^M ifandonlyif S (wW*i 

V»2/ V2/2 / v^2/ Vy2 



The operation of obtaining B from A has been encountered in several contexts. Tucker [15] 
considers an equivalence relation among rectangular matrices, which is implicitly determined 
by a nonsingular (not necessarily principal) submatrix and is defined as follows: two m by n 
matrices A and C are combinatorially equivalent if there is a one-to-one correspondence between 
the sets of ordered pairs {(x, y)\y = Ax} and {(u, v ) | v = Cu}, given via a permutation matrix 
P of order m + n by (x T ,y T ) T = P (u T ,v T ) T . It is shown in [15] that C is combinatorially 
equivalent to A in (|l.l| ) if and only if C is, up to independent permutations of its rows and 
columns, equal to B in (|l.2| ) with the signs of the off-diagonal blocks reversed. The matrix C 



is referred to as a pivotal transform of A. When the equivalence relation is determined by a 



principal submatrix, Tucker [16] refers to C as a principal pivotal transform of A and asserts 
that if A has positive principal minors (that is, if A is a P-matrix), then so does every principal 
pivotal transform of A (see Theorem |5.3| ). 



In the sequel we will adopt the more commonly used term of 'principal pivot transform'. 

Tucker's motivation for introducing combinatorial equivalence and studying principal pivot 
transforms is rooted in an effort to generalize Dantzig's simplex method from ordered to gen- 
eral fields. In turn, the domain-range relation between A and B observed by Tucker is later 
used by Cottle and Dantzig |5[ as an important feature of their "principal pivoting algorithm" 
for the linear complementarity problem when the coefficient matrix is a real P-matrix. In that 
algorithm, principal pivot transforms are used to exchange the role of basic and nonbasic vari- 
ables of the problem and the fact that principal pivot transformations preserve P-matrices is 
applied effectively. Principal pivot transforms have since found similar uses in the context of 
mathematical programming (see e.g., Pang [|k|]). 



The relation between A and B above prompted Stewart and Stewart [14] to refer to B as the 
exchange of A (exc(A)). The authors use exchanges in order to generate S-orthogonal matrices 
from hyperbolic Householder transformations, and then apply them to solve the mixed Cholesky 



updating/downdating problem. In [14| it is also noted that this method of construction of S- 
orthogonal matrices is a folk result in circuit theory and a reference to Belovitch [||] is made for 
a special case. 
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In Johnson and Tsatsomeros [11], a fundamental matrix factorization of the principal pivot 



transform turns up in a discussion of row- interval nonsingularity and the relation to P-matrices. 



We review this factorization in Lemma 3.4 and take the opportunity to provide a proof valid 



for complex matrices of a result claimed in |TT[ (see Remark 5.4). In a related vein, Eisner and 



Szulc H introduce a generalization of P-matrices to block P-matrices and show that a certain 
class of block P-matrices is left invariant under principal pivot transformations. 

The principal pivot transform also appears under the term gyration in Duffin, Hazony, and 
Morrison 0, and is mentioned in a survey of Schur complements by Cottle g. 

The above varied interest for principal pivot transforms motivates us here to survey and further 
study their general properties. We will discuss the determinants, the eigenvalues and other basic 
characteristics of principal pivot transforms relative to arbitrary principal submatrices. The 
relation and parallelism of the principal pivot transformation to inversion will also be considered, 
as well as a potential application to iterative techniques for solving linear systems (see sections 3 
and 4). We will also discuss matrix classes left invariant under principal pivot transformations, 
including the aforementioned P-matrices and S-orthogonal matrices (see section 5). 



2 Notation and preliminaries 

Let n be a positive integer and A G M n (C). The i-th entry of a vector x is denoted by x(i). In 
the remainder the following notation is also used: 

• ( n ) = {1)2,..., n}. For any a C (n), the cardinality of a is denoted by |a| and a c = a\(n). 

• A[a, 0\ is the submatrix of A whose rows and columns are indexed by a, (3 C (n), respec- 
tively; the elements of a, (3 are assumed to be in ascending order. When a row or column 
index set is empty, the corresponding submatrix is considered vacuous and by convention 
has determinant equal to 1. 

• A[a] = A[a,a], A(a,(3] = A[a c ,(3]; analogously we define A[a,(3), A(a,/3) and A(a). 

• A/A[a] is the Schur complement of an invertible principal submatrix A[a] in A, namely, 
A/A[at] = A(a,a) — A(a, a](A[a, a])~ 1 A[a, a). It is well known that det(^4/A[a]) = 
det A I det (A [a]). 

• a{A) is the spectrum and p(A) the spectral radius of A. 

• diag(di, . . . , d n ) is the diagonal matrix in M n (C) with diagonal entries d%, . . . , d n . 



Definition 2.1 Given a C (n) and provided that A[a] is invertible, we define the principal 
pivot transform of A € M n ((D) relative to a as the matrix ppt (^4, a) obtained from A by 
replacing 
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A[ct] by A[a] _1 , A[a,a) by -A[a]~ l A[a,a), 
A(a, a] by A(a, a]^4[a] _1 , and A{a) by A/A[a]. 

By convention, if a = 0, then ppt (A, a) = A. 

The principal pivot transform is related but distinct from the following block representation of 



the inverse (obtained by combining formulas in [|3| and |18[; see also [10, section 0.7.3]): Given 
an invertible A £ M n ((D) and a C (n) such that A [en] and j4(a) are invertible, A -1 is obtained 
from A by replacing 

A[a] by (A/A(a))~ l , A[a,a) by -A[a]"M[a, a){A/A[a})- 1 , 
A(a,a] by (A/A[a]) _1 A(a, a]A[a] _1 , and A(a) by (A/A[a]) _1 . 

In our subsequent discussion, we will also use an easy to verify determinantal formula for A + D, 
where D = diag(di, . . . , d n ), namely, 

det(A + D) = Y[didetA[a}. (2.1) 

aC{n) i^a 

3 Basic properties of principal pivot transforms 

We begin with a formal statement of the basic domain-range exchange property of ppt (A, a), 
and include a proof sketch for the sake of completeness. 

Theorem 3.1 Let A £ M n ((D) and a C (n) so that A[a] is invertible. Given a pair of vectors 
x,y £(D n , define u,v £(D n by u[a] = y[a], u(a) = x(a), v[a] = x[a], v(a) = y{a). Then 
B = ppt (A, a) is the unique matrix with the property that for every such x,y, y = Ax if and 
only if Bu = v. Moreover, ppt (B,a) = A. 

Proof. Consider the permutation matrix P for which 

Px=[ X[a] ) and PAP T -( A[a] A[aM 



.x(a) J \A(a,a] A(a) 

By the construction outlined in Definition ETA and on letting B = ppt (^4, a), we have 



PBP 



■T 



Ala}^ 1 -A[a]~ l A[a.,a) 
A(a,a]A[a]~ 1 A/A[a] 



Then, with u and v as prescribed, it can be easily verified that PAP T (Px) = Py if and only 
if Pv = PBP T (Pu), or equivalently, Ax = y if and only if Bu = v. To show uniqueness, 
suppose that B'u = v if and only if Ax = y. Then (B — B')u = for all u such that 
u[a] = y[a] = A[a]x[a] + A[a,a)x(a) and u(a) = x(a). As A[a] is invertible and x is chosen 
freely, it follows that (B - B')u = for all u £© n , that is B = B' . To see that ppt (B, a) = A, 
notice that ppt (B, a)x = Ax for all x € C n . a 
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It is interesting to note in the next theorem that in certain cases, consecutive principal pivot 
transforms result into the inverse of a matrix. 

Theorem 3.2 Let A G M n ((D) and suppose that there exists a partition of (n) into subsets 014, 
i = 1, 2, . . . , k so that the sequence of matrices 

A = A, Ai = ppt(A i ^ 1 ,a i ), i = l,2, ...,k 

is well defined (i.e., the matrices Ai—i[ai\ are invertible). Then A is invertible and A" 1 = A^. 



Proof. By Theorem 3.1 applied to each of the A^ in sequence, and since the on are mutually 
disjoint and their union is (n), we have that Ax = y if and only if A^y = x for all x,y G(D n . It 
follows that A is invertible and by uniqueness of the inverse that A^ = A^ 1 . m 

Remark 3.3 In IE] it is observed that A -1 G M n ((D) can be found with a sequence of at most 
n principal pivot transforms (and by interchanging rows or columns if needed). Adopting the 
definition of a flop as the time required to execute x = x + t*x, we can compare such a method 
of inversion of A G M n (lR) with solving the n linear systems Ax = ej via the LU factorization of 
A. The latter method of inversion of A entails n 3 /3 + n • n 2 /2 = 5n 3 /6 flops. Suppose now that 
the partition 014 = {i}, i = 1, 2, . . . , n of (n) yields, as in Theorem |3.2| , the inverse of A. In the 
process, it is firstly required to compute A\ = ppt(A,a\), which entails 2(n — 1) divisions (for 
the off-diagonal blocks), (n — l) 2 multiplications for the computation of the Schur complement 
(which is a rank one update of A(a)), and 1 division for the calculation of 1/an. The total 
is therefore n 2 flops for the calculation of A\. Thus to find the inverse by calculating A n , the 
required flop count is 

n 2 + (n-l) 2 + ... + 2 2 = n(n + 1)(2re + 1) -l. 

6 

It follows that there is an economization of (3n 3 — 3n 2 — n + 6)/6 flops over inversion via LU 
factorization that can be realized e.g., when the inverse of a P-matrix is sought (cf. ^.2[). 



To study further the basic properties of ppt (A, a), we continue with a useful observation that 
appears implicitly in the proof of |I5|, Theorem 4] and in [|ll]] . 

Lemma 3.4 Let A G M n ((D) and a C (n) so that A[a] is invertible. Let T\ be the matrix 
obtained from the identity by setting the diagonal entries indexed by a equal to 0. Let T2 = I—T\ 
and consider the matrices C\ = T2 + T\A, C2 = T\ + T2A. Then ppt (A, a) = CiC^ 1 ■ 

Proof. Without loss of generality, we can assume that a = (k) (otherwise we can apply our 
argument to a permutation similarity of A) . Observe then that 

<*-(,.' , ° ) A ^ 

\A[a, a\ A{a) J \ I 
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and thus 



C\C 2 



-i 



/ 

A(ot, a] A(a) 

A(a,a]A[otf l 



Ala]' 1 -A[a]~ 1 A[a,a) 
/ 



-A[a]~ l A[a,a) 
A/A[a] 



ppt (A, a). 



Definition 3.5 Referring to the matrices of Lemma 3.4, we call ppt (A, a) = C\C 2 1 the basic 
factorization of ppt (A, a). 



In connection with a remark added in proof in [15], we have the following result that sheds more 
light on the combinatorial relationship between a matrix and its principal pivot transforms. 



Theorem 3.6 Let A £ M n ((D) and a C (n) so that A[a] is invertible. Let B = ppt {A, a) and 
I be the identity in M n ((D). Then there exists a permutation matrix P £ M2 n ((D) such that 



i-B I r [ l A ) 



(3.1) 



Moreover, if T\ and T 2 are as in Lemma ^L£, then 

(T x T 2 



\T 2 Tx)' 



P 



Conversely, if B = ppt (A, a) for some a C (n), then (|3.l| ) holds for an appropriately defined 
permutation matrix P £ M 2n (<E). 



Proof. In the notation of Lemma 3.4, we have that B = ppt (A, a) if and only if 



(-B I) 



Co 

Ci 



0. 



The claims of the theorem follow by substituting C\ = T 2 + T\A and C 2 = T\ + T 2 A. That P 
as above is a permutation matrix follows from the fact that T\ +T 2 = I. _ 



Example 3.7 To illustrate the definitions and observations so far, let a = {1,3} so that 



A 



fl 2 1' 

1 1 ] and B = ppt {A, a) 
\2 8 1, 




Notice the exchange taking place relative to the index set a in the equations 
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The basic factorization of B is C\C 2 1 , where 



Ci 




and Co 



1 


2 


1 





1 





2 


8 


1 



Also if p = {2}, then 



ppt (B,f3) 



{ 


.2 


1.2 


-.2 




-.2 


-.2 


.2 


\ 


1.2 


-.8 


-.2 



A' 



Theorem 3.8 Let A £ M n ((D) and a C (n) so i/iaf .A [a] is invertible. Then 

... , . , . detA(a) 
(%) det(ppt (A, a)) = , and 

det A[a\ 

(ii) if in addition A(a) is invertible, ppt (A, a) -1 = ppt (A, a°). 

Proof. Let CiC^ 1 be the basic factorization of ppt (A, a). The conclusions follow, respectively, 
from Lemma |3A| and by directly verifying that ppt (A, a)" 1 = C^C^" 1 = ppt (A, a c ). B 

Note that invertibility of ppt (A, a) does not necessarily imply invertibility of A (Meenakshi 
fli~2|| ). A simple counterexample is provided by 

A=(] and ppt(A{l}) = (J 



4 Eigenvalues of principal pivot transforms 

We continue with what to our knowledge are new observations on the eigenvalues of principal 
pivot transforms. 

Theorem 4.1 Let A £ M n ((D) and a C (n) so i/iai A[a] is invertible. Let CiC^ 1 be the basic 
factorization of ppt (A, a). Then the following are equivalent: 

(i) A G o-(ppt (A, a)) 

(ii) A is a finite eigenvalue of the matrix pencil C\ — AC2. 

When, in addition, A 7^ 0, then the following condition is also equivalent to (i) and (ii): 

(Hi) A — D is singular, where D = diag(di, . . . , d n ) with di = A -1 if i £ a and 
di = A otherwise. 
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Proof. The equivalence of (i) and (ii) follows from Lemma 3.4 and the fact that A is a finite 
eigenvalue of the matrix pencil C\ — XC2 if and only if A is an eigenvalue of CiC^ 1 . For the 
equivalence of (ii) and (iii) when A 7^ 0, observe that up to a permutation similarity of A, 



Ci — AC2 



I[a]-AA[a] -XA[a,a) 
A(a,a] A(a) — XI (a) 



where I is the identity matrix in M n ((D). Thus, multiplying the leading \a\ rows of C\ — XC2 
by —A -1 , we obtain that (ii) holds if and only if 



A" 1 I[a] 
A 1(a) 



is singular. 



It is worth noting the parallelism in viewing a principal pivot transform as 'partial inversion' 
with the fact that its nonzero eigenvalues are the zeros of det(A — D) as in (iii) of the above 
theorem. A more precise account of det(A — D) as a function of A and of its relation to the 
spectrum of the principal pivot transform is given next. Note that unless a = 0, det(^4 — D) is 
not a polynomial in A. 

Proposition 4.2 Let A G M n ((D) and a C (n) such that A[a] and A(a) are invertible. Let X 
be an indeterminate, and let D = diag(c?i, . . . ,d n ) with di = A -1 if i G a and d, L = X otherwise. 
Then 

w det(,4 - D) 
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(A) = (-l)l QC lA |a| 



detvlfal 



is the characteristic polynomial o/ppt (A, a). Moreover, the coefficients of g(X) can be expressed 
as real linear combinations of the principal minors of A. 

Proof. Since A[a] and A(a) are invertible, we respectively have that B = ppt (^4, a) is well 
defined and, by Theorem |3.8| , nonsingular. It then follows from Theorem |4.1| (iii) that A is 
an eigenvalue of B if and only if det(^4 — D) = 0, where D is as described above. Since D is 



diagonal, by (^1) we obtain 



det(A-D)= 53 (-1)1^1 JJdjdetA^. (4.2) 

/3C(n> m 

Since (3 C = {(3 C n a c ) U {(3 C n a) and (fj c n a c ) n {(3 C n a) = 0, we have 

Y[ dl = 5] Al" w H^I. (4.3) 

Also notice that \a\ > \/3 c Da\> \(3 C n a\ - \(3 C n a c \, i.e., 

\{3 c na c \ - \(3 c Da\ > -\a\, (4.4) 

and that 

\(3 C n a c \ - \(5 C n a\ < \(3 C n a c \ < \a c \. (4.5) 
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Equalities hold in ( [Q| ) and ( |4.5| ) if and only if /3 C = a and /3 = a, respectively. Thus, multiplying 
the equation in Q4.2|) by Al Q and using (4.3)-(4.5), we obtain that A is an eigenvalue of B if and 



only if A is a (nonzero) root of the polynomial 



Al Q l det (A -D)=Y, (-l^lA^+l^^l-l^l det A[f3}. (4.6) 

/9C(n> 



The term of highest degree in (|4.6| ) appears when (3 = a and equals (—1)'° A n detj4[a]. The 
constant term in fl4.6| ) appears when /? = a c and equals (—1)'°' det A(a). Thus, by Theorem 
|3.8| (i), g(X) as in the statement of the theorem is indeed the characteristic polynomial of B and 
its coefficients are real linear combinations of the principal minors of A as seen by 



Note that under the assumptions (and as a consequence) of the above proposition, if A G M ra ((D) 
has real principal minors, then the spectrum of ppt (A, a) is closed under complex conjugation. 



Corollary 4.3 Let A G M n ((D) and a C (n) so i/iai A [a] is invertible. Then 1 € <r(ppt (A, a)) 
(resp., —1 £ o" (ppt (A, a)) ) if and only if 1 G o"(A) (resp., —1 G c(A),). ^4Zso ppt (A, a) is 
singular if and only if A(a) is singular. 



Proof. The results on the ±1 eigenvalues follow from Proposition 4.2. The singularity condition 



for ppt (A, a) follows either from Theorem 3.8 (i) or from Theorem |4.1| (ii). 



We continue with an application to iterative techniques for solving a linear system Ax = b, 
where A G M„(C) is invertible. Such iterative techniques are obtained by expressing the 
unique solution x as a fixed point of a matrix equation x = Tx + c for an appropriate matrix 
T. In fact, based on a splitting of A into A = M — N and assuming that M is invertible, we 
take T = M _1 iV and c = M~ x b. Then the sequence {xk}^ generated by Xk = Txk-i + c for 



arbitrary xq converges to the solution x if and only if p(T) < 1 (see e.g., Varga [17]). The 
Jacobi method is obtained when M = diag(an, . . . , a nn ) and N = M — A. In many instances, 
certain splittings lead to divergent sequences. This may be overcome by considering a principal 
pivot transform T of T and an equation x = Tx + d equivalent to x = Tx + c, as suggested by 
the following result and illustrated by the subsequent example. 



Proposition 4.4 Let T G M n ((D) and x,c G(D n . Let a C (n) so that T[a] is invertible. 
Consider the vector u defined by 




c(i) if i G a 
otherwise. 



Then x = Tx + c if and only if x = Tx + d, where d = c — (L + T)u. 



Proof. Let T, T, x, c, u and d as prescribed. Observe that by Theorem 3J, Tx = x — c is 
equivalent to T(x — u) = x — (c — u), which in turn is equivalent to x = Tx — Tu + (c — u), that 
is, x = Tx + d. m 
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Example 4.5 Consider the matrix A and the corresponding Jacobi iteration matrix T given 
by 

/ 3/2 1/4* 
and T= 3/2 
Vl/2 1/2 

That is, as p(T) > 1, the Jacobi iteration 
Txk-i + c fails to converge to the solution of a system Ax = b. However, if we consider 





I 1 


-3/2 


-l/4\ 


A = 


-3/2 


1 


Tj 




V-l/2 


-1/2 




We find that a(T) = 


{2.1419, 


-.6419, 


-1.5}. 



f = p P t(r,{i,2}) 



/ 2/3 

2/3 
Vl/3 1/3 




then <r(T) = {—1/4,0,2/3} and thus p(T) = 2/3 < 1. It follows that the iteration = 
Txk-i + d with d as in Proposition ^4|, converges to the solution of Ax = b. In passing we 
mention that T above satisfies the assumptions of the Stein-Rosenberg theorem in [17| and 
hence the Gauss-Seidel iteration for A also fails to converge to the solution of the system. 



5 Principal pivot transforms of special matrices 

One of the main matrix classes discussed in association with principal pivot transforms is the 
class of P-matrices, that is, matrices in M n ((D) all of whose principal minors are positive. 
Tucker |uj asserts that principal pivot transformations preserve the class of P-matrices. In the 
case of real P-matrices, a simple proof of this assertion can indeed be based on Theorem |3.1| 
and on the following characteristic property of real P-matrices (see Fiedler g, Theorem 5.22]): 
A E M n (H) is a P-matrix if and only if for every nonzero x E IR n , x and Ax have at least one 
pair of corresponding entries whose product is positive. Here we present a proof of the assertion 



in [16 1 for the general case of complex P-matrices, based on the following well known result. 



Lemma 5.1 Let A E M n ((D) be a P-matrix and a C (n). Then A/A[a] is a P-matrix. 

Proof. Assuming that A is a P-matrix and by considering the block representation of A^ 1 
mentioned in section 2, it is enough to show that A' 1 is also a P-matrix. Indeed, since A[a] is 
invertible for all a C (n), each principal submatrix of A -1 is of the form (^/^[a]) -1 for some 
a E (n) and thus its determinant is det A [a]/ det A > 0. B 

Theorem 5.2 Let A E M n ((D) be a P-matrix and a C (n). Then ppt (A, a) is a P-matrix. 

Proof. Let A be a P-matrix and consider first the case where a is a singleton; without loss 
of generality assume that a = {1}. Let B = ppt (A, a) = (&«)■ By definition, the principal 
submatrices of B that do not include entries from the first row of B coincide with the principal 



submatrices of A/ A [a] and thus, by Lemma 5.1, have positive determinants. The principal 



submatrices of B that include entries from the first row of B are equal to the corresponding 
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principal submatrices of the matrix B' obtained from B using bu = (A[ot\) 1 > as the pivot 
and eliminating the nonzero entries below it. Notice that 



B' 



1 0\ / b\\ — b\\A\a, a) \ /fen —bnA[a,a) 
-A(a,a] I J \A(a,a]b n A/A[a] J ~ \ A(a) 



That is, B' is itself a P-matrix, as it is block upper triangular with the diagonal blocks being 
P-matrices. It follows that all the principal minors of B are positive and thus B is a P-matrix. 
Next, consider the case a = . . . ,ik} Q (n) with k > 1. By the proof completed so far, the 
sequence of matrices 

A = A, Aj = ppt (Aj-i, {ij}), j = 1, 2, . . . , k 
is well defined and comprises P-matrices. Moreover, from the uniqueness of B = ppt (A, a) 



shown in Theorem 3.1, it follows that A^ = ppt (A, a) = B and thus B is a P-matrix. 



The next theorem summarizes our discussion of principal pivot transforms of P-matrices and 
follows readily from the above result. 

Theorem 5.3 Let A € M n ((D). Then the following are equivalent: 

(i) A is a P-matrix. 

(ii) there exists a C (re) such that ppt (A, a) is a P-matrix. 
(Hi) for all a C (n), ppt (A, a) is a P-matrix. 



Remark 5.4 Theorem 5.3 is stated in similar terms in [11, Theorem 4.4]. However, the proof 



provided in [11|, unless modified, is valid only when A is a real matrix. 



We continue with a few words on some other matrix classes that are invariant under principal 
pivot transformations. One such class is the S-matrices or semipositive matrices, consisting 
of matrices A E M n (IR) such that Ax > for some x > (inequalities here are entrywise.) 
Clearly, by Theorem |3.1|, a principal pivot transform of an S-matrix is an S-matrix. 



Next, recall that A € M n (JR) is called a Z-matrix if its off-diagonal entries are all nonpositive. Of 
course, principal pivot transformations do not, in general, preserve Z-matrices. In particular, 
they do not preserve M- matrices (i.e., Z-matrices that are also P-matrices; see Berman and 
Plemmons [0]). However, principal pivot transformations do preserve a class that generalizes 



M-matrices, which is introduced in [13|. The matrix A € M n (IR) is called a hidden Z-matrix 
provided there exist Z-matrices X, Y such that 



AX = Y and r T X + s T Y > 



for some vectors r, s > 0. As it is shown in 13], principal pivot transformations preserve the 
intersection of the classes of hidden Z-matrices and P-matrices. For example, the principal 
pivot transform of an M-matrix is a hidden Z-matrix and a P-matrix. 
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We now return to the S-orthogonal matrices mentioned in the introduction. The matrix Q G 
M n (]R) is called S-orthogonal if there exists a signature matrix S G M n (H) (that is, a diagonal 
matrix S whose diagonal entries are ±1) such that Q T SQ = S. When S = I, then an S- 
orthogonal matrix is simply an orthogonal matrix. In []l4|] it is formally shown that S-orthogonal 
matrices can be constructed for any prescribed signature matrix S in the following way. Suppose 
that S = diag(si, . . . , s n ) and that Sj = 1 for all i G a C (n) and = —1 for all i G a c . Let 
i? G M n (H) be an orthogonal matrix such that R[a] is invertible. Then Q = ppt (R,a) exists 
and is S-orthogonal. 

As is the case with Schur complements, the notion of a principal pivot transform can be extended 
to the case of non-invertible principal submatrices by considering generalized inverses. Some 
work in this direction is presented in fL^f , where it also shown that under certain assumptions, 
the principal pivot transform of an EP-matrix is an EP-matrix. ( Recall that A G M n ((D) is an 
EP-matrix if Nul(A) =Nul(A*).) 

6 Some questions 

We conclude with a couple of questions about principal pivot transforms, hoping to motivate 
their further theoretical development and to promote their applicability. 

It has been shown in Coxson || that the important problem of testing for P-matrices is co-NP- 
complete. In view of Theorem |5.3| , we are led to ask: Is there a computationally advantageous 
utilization of principal pivot transforms to check whether a given matrix is a P-matrix or not? 

As we saw in section 4, principal pivot transformations in certain instances can map the eigen- 
values to desired regions, e.g., the open unit disk. When and how can we choose a so that the 
eigenvalues ofppt(A,a) lie in given regions of the complex plane? 
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